Absence of pre-classical solutions in Bianchi I loop quantum cosmology 
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Loop quantum cosmology, the symmetry reduction of quantum geometry for the study of various 
cosmological situations, leads to a difference equation for its quantum evolution equation. To ensure 
that solutions of this equation act in the expected classical manner far from singularities, additional 
restrictions are imposed on the solution. In this paper, we consider the Bianchi I model, both 
the vacuum case and the addition of a cosmological constant, and show using generating function 
techniques that only the zero solution satisfies these constraints. This implies either that there 
are technical difficulties with the current method of quantizing the evolution equation, or else loop 
quantum gravity imposes strong restrictions on the physically allowed solutions. 
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When a new method of looking at physical situations 
arises, it is frequently too complicated to look at the gen- 
eral case, and thus it is necessary to look at simpler ver- 
sions as test models. For example, in quantum mechanics 
one considers the hydrogen atom before examining more 
involved systems. These reduced models not only are 
easier to solve, but frequently provide information about 
what to expect from the more generic situation. This is 
the case with loop quantum cosmology 0], a symmetry 
reduced version of the full theory of quantum geome- 
try 0] (also known as loop quantum gravity), which is 
intended to provide a framework joining together general 
relativity and quantum mechanics. 

Much like the Bohr quantization condition for the elec- 
trons of the atom, space-time itself is no longer continu- 
ous in quantum geometry, but instead is discrete. Since 
this means the states representing geometry are also dis- 
crete, the quantum evolution equation becomes a differ- 
ence equation. In general, however, solutions to a re- 
cursion relation will oscillate in sign as one moves in a 
particular direction in the parameter space. Ordinarily, 
unphysical states like these would be eliminated by the 
imposition of an inner product, but this is unavailable at 
the moment in loop quantum cosmology. To get around 
this, a notion of pre-classicality |3( is used to pick out 
those quantum solutions to match a semi-classical state 
far away from the singularity. This condition will elimi- 
nate all states which vary greatly at Planck-length scales. 
We will see here how the pre-classicality of a solution can 
be imposed as a global property. 

Not only does loop quantum cosmology allow us to 
study simple situations of physical interest - for exam- 
ple, the quantized version of the isotropic Friedman- 
Robertson- Walker (FRW) space-times that serves as the 
basis of modern cosmology - so that physical predictions 
can be made, but it can give information about the full 
theory. In particular, it is closer in spirit to full quantum 



geometry than the minisuperspace models considered in 
standard quantum cosmology would be to a quantized 
general relativity. For the latter, the metric was first re- 
duced to a specific model, then the equations of motion 
were quantized. Loop quantum cosmology takes the op- 
posite path: starting with the kinematic space of quan- 
tum states from the full model, those states that obey 
the particular symmetry are chosen and then the Hamil- 
tonian constraint on these states is quantized. 

So far, most work on loop quantum cosmology has 
been done on isotropic models, looking at both quantum 
and semi-classical aspects. This has included a resolu- 
tion of the classical singularity [J], cosmological appli- 
cations such as inflation [5| and oscillatory universes |fj, 
and semi-classical approximations of the quantum Hamil- 
tonian constraint jjj. However, some anisotropic models 
have also been studied, including Bianchi IX 8] and the 
Bianchi I model |2( when a local rotational symmetry 
(LRS) is assumed. This allows us to see how loop quan- 
tum cosmology fares in more general settings, while still 
being simple enough to get exact solutions to the quan- 
tum Hamiltonian constraint. In particular, for Bianchi 
I LRS, generating function techniques 01 were used to 
find possible solutions to this constraint, a partial differ- 
ence equation in two parameters. By placing boundary 
conditions on the differential equation associated with 
this recursion relation, one could find all pre-classical so- 
lutions to the original quantum evolution constraint. In 
this paper, we examine the full Bianchi I model and carry 
out the same program (this model has also been looked at 
using related methods [111). Yet, as we will see, the result 



we reach will be a negative one 
solutions arc the null solutions. 



the only pre-classical 



The Bianchi I model is one of the simplest anisotropic 
models, with a metric of the form ds 2 = —dt 2 + 
J2i a i(t)dx 2 {i — 1,2,3), where the functions a,k(t) are 
the scale factors in the respective spatial direction. The 
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flat (A; = 0) FRW metric used in cosmology is the spe- 
cial case where the three scale functions are identical. 
In the connection formalism [l^ . there are nine invari- 
ant degrees of freedom, which are reduced to six after 
the gauge freedom is removed. When restricted to di- 
agonal metrics, we have that the connection 1-form is 
given by A z a = cr K \A l K u)^ , where (because of the diag- 
onalization) are a fixed set of left-invariant 1-forms, 
and A l K 6 SO(3) is a rotation matrix. The momenta 
conjugate to the coefficients C(k) ar e given by the compo- 
nents of an invariant densitized triad Ef = p( K > Af X 1 ^-, 
where the X^ are the vector fields dual to ui^ . These 
triad components are related to the scale factors cii of a 
Bianchi I metric by p\ = \a2a3\sgn(ai), and so on for the 
other two momenta. The symplectic structure given by 
{ci,p J } = KjSj, where 7 is the Barbero-Immirzi param- 
eter of loop quantum gravity and n = 8irG. 

We now briefly review the quantization of the Bianchi 
I model 0| ■ Because of the diagonalization, the basis 
states will be three copies of those used in the isotropic 
model. Starting with the eigenstates of the triad 

ra) = — = 

V2sin(f) 

where m is a real number, a gauge-invariant state |s) can 
be expanded as a series 
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|mi) <g> |m 2 ) ® \m 3 ), 



mi ,7772 ,m3 



where the coefficients s 



mi ,7772 ,7773 



satisfy the conditions 



,7772 ,7773 



-mi, — 7772,7773 



3 mi , — 7772 , — 7773 ' 



-mi 5 m 2 1 — m 3 



These latter conditions come from the remaining freedom 
to change the sign of two of the triad vectors simultane- 
ously; it allows us to choose one of the three parameters 
to run from —00 to 00, and restrict the other two to non- 
negative values. For example, if we let 7773 S (— 00,00), 
then TO3 = would correspond to the classical singu- 
larity. In this work, however, we will look only at the 
subset where all three parameters are non-negative; as 
we will comment below, the symmetry of the Hamilto- 
nian constraint would allow one to build up a solution 
for all possible values of nik- 

When this classical structure is carried over into quan- 
tum operators, we get a partial difference equation for 
the Hamiltonian constraint in terms of the modified co- 
efficients t mi , m2 ,m 3 = V2mi,2m2,2m 3 S2mi,2m 2 ,2m, 3 , where 

Kni,m 2 ,m 3 = {^^ 2 p) 3/2 Vl m i m 2 m 3 1 are the eigenvalues 
of the volume operator V acting on the state |mi) <g) 
|m 2 ) <8> \rri3}, and tp is the Planck length. Because of 
their definition, with ( mijm!iraj including a volume fac- 
tor, we have that, e.g. io. 7712,7773 = 0, and similarly for the 
other boundaries. We can use the sequence t since the 



original values so, 7779.7773? ^mi,o,m 3 and s mi rn2 o fall out of 
the recursion relation and thus have no effect on the final 
solution. Also, the parameters are all real numbers, 
but the difference equation relates only those separated 
by a jump of one step. Here, we will examine those se- 
quences that pass through the classical singularity; the 
rest of the sequence can be built up by continuity. Be- 
cause the intermediate values will be interpolated, there 
are no problems scaling the parameters nik when map- 
ping between the s and t sequences. If we define the 
difference operator <5i as 



t 



7711+1,7712)7713 ^7/71 — 1,7772 , 7773 5 



and similarly for 62 and 63, the Hamiltonian constraint 
in Bianchi I for the particular case of a cosmological con- 
stant A is given by the difference equation 



[d(mi)S2S 3 + d(m,2)5iS 3 



mi , rn 2 , m 3 



d(m 3 )S 1 6 2 
= 0, 



(1) 



where 



d(n) 



1 

2r> 



1 

2m 



Partial difference equations are more difficult to solve 
than recursion relations of one parameter, so next we 
discuss how to find solutions. 



In order to find a sequence t 



7771 ,7772 , 7773 



that satisfies the 
recursion relation Q , we will use the method of gener- 
ating functions, briefly described here. For more details 
using this technique in loop quantum cosmology, see 
We will define a function of three variables such that 



F(x,y,z) 



$3 u "< 1 

HI | — 7772—0 777 3 =0 



x mi ^7772 .,7773 ^ 



and find a partial differential equation for F(x, y, z) that 
is equivalent to the recursion relation for t mi ,m2,m 3 - It 
is important to realize that finding a solution F gives 
the entire sequence t - there is no need to choose a par- 
ticular parameter as a "time" to evolve in, avoiding 
such issues of interpretation. Notice that we only include 
positive powers of the variables x, y and z in the series, 
whereas one of the parameters, say, 7773, can assume neg- 
ative values, passing through the singularity 7773 = 0. 
Because of the symmetry of the recursion relation 
under 7773 — > —7773 from the fact that d(—m 3 ) = —d(m 3 ), 
we can match any two of our solutions at this singularity. 

The question now is how to pass between a recursion 
relation for a sequence, and a condition on the associated 
generating function. Suppose we have a one-parameter 
sequence Tk and corresponding generating function R(x). 
Multiplication by a power x d will result in shifting the 
sequence — > rk+d, while multiplication by a function 
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of the parameter k, e.g. kr^, is related to differentia- 
tion by x [ljj. By using ideas of this type, we can find 
a differential equation for F(x, y, z) which has the same 
information as the recursion relation ft). However, the 
generic solution found by this method will not be pre- 
classical, but instead will oscillate in sign as one moving 
towards increasing to&; we must invoke additional con- 
ditions on the function F so that we only get physical 
solutions. 

To do this requires some thought about which se- 
quences will avoid sign oscillation, ft turns out that this 
depends crucially on the pole structure of the function F. 
As an example, a simple pole at x = f will lead to a Tay- 
lor series expansion whose coefficients become constant 
asymptotically, since (x— f) _1 = 1 + x + x 2 +x 3 + ■ ■ ■ . On 
the other hand, we wish to avoid simple poles at x = — 1, 
since the coefficients of the Taylor series alternate sign, 



with (x + ly 1 = 1 - x + x 2 



So having a 



pole at x — 1 will result in a sequence whose coefficients 
are asymptotically constant. However, to avoid the se- 
quence growing without bound, this must be a simple 
pole; this behavior can be seen in the expansion of the 

double pole (x - f )~ 2 = 1 + 2x + 3x 2 + 4x 3 H . Again, 

because we are able to look at the global features of the 
sequence i mijm2 , m3 by fixing the properties of the func- 
tion F, we do not have to postulate a particular wave 
form that the sequence must "evolve" to at large values 
of a parameter m/j. Instead, we can stipulate that the se- 
quence must asymptotically converge on a solution which 
has the desired pre-classical properties; the model itself 
tells us what kinds of asymptotics are possible. 

With these ideas, we start with the Hamiltonian con- 
straint ft): however, because the function d(n) is not 
polynomial, we cannot easily change functions of m& into 
differential operators. So we make an approximation to 
allow us to find a differential equation for the generating 
function; using the fact that d(n) ~ l/2n, we can write 
the recursion relation (|TJ as 



(■mim 2 5iS 2 + mim 3 SiS 3 + m 2 m 3 525z 



2K7 3 £A)i 



mi ,7712 ,™-3 



0. 



(2) 



Since d(0) = 0, this equation works only for 
m i,m2,m3 > 1. Because we are making an approxima- 
tion for the function d(ri), we should have terms on the 
right-hand side of this recursion relation with factors of 
the difference [d(rrik) — \ j2m^\. However, the size of the 
error if these terms are neglected will be small compared 
to the sequence itself [ij. 

By using the mapping between the recursion relation 
and partial differential equations discussed above, we can 
find an equation for a three- variable function F(x,y,z), 
whose solutions correspond to sequences t 



mi ,m2 ,m3 



that 

solve the approximate recursion relation J3) for the 
Bianchi I model. This is given by 



d 2 



(\-x 2 ){\-y 2 )F 



dxdy 
d 2 



dydz 



(1 



xyz 



(I - x 2 )(l - z 2 )F 



z 2 )F 



xyz 



dxdz |_ xyz 

1 dxdydz 



(3) 



With this equation, it is natural to define a new function 
H(x, y, z), such that 



F(x,y,z) 



xyzH(x,y,z) 



(1 - x 2 )(l - y 2 )(l ~ z 2 ) 



Not only does this make the differential equation © 
for the generating function more tractable, but it also 
gives us a way to ensure that the function has the 
proper asymptotic behavior. By assuming that the func- 
tion H(x,y,z) has boundary conditions H(—l,y,z) = 
H(x,—l,z) = H(x,y,—1) = 0, the original function 
F(x,y, z) can have at most a simple pole at x = — 1, etc. 
The requirement for a pre-classical sequence f mi , m2 . m3 is 
manifest as a set of boundary conditions for the derived 
function H. When we substitute for F in terms of H, 
and multiply by (f — x 2 )(l — y 2 )(l — z 2 ), this gives us the 
equation 



(l-x 2 )(l-y 2 )- 



+ (1 - x z )(l - z z ) 



dxdy 

-2 KJ %A(l-x 2 )(l-y 2 )(l-z 2 ) 



dxdz 
d 3 



dxdydz 



x 2 y 2 z 2 H 
(f - x 2 )(l - y 2 )(l - z 2 ) _ 



(4) 



with boundary conditions mentioned above. We show H(x,y,z) = 0. 
next that, regardless of the value of the cosmological 

constant A, this equation only has the trivial solution Because of the boundary conditions that are necessary 

to assure pre-classicality, any solution H(x 7 y,z) will be 
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the sum of terms in the form (x + \) l (y + iy(z + l) fe - 
We will use this fact to show that the case A = has 
only trivial pre-classical solutions. First, we define new 
variables u = x + l,v = y + l, and w = z + 1, and look 
at the lowest order terms in u,v and w, namely 



H(u, v, w) 



i+j+k—d 



Cij,kU l v j w k 



with the parameter d defining the order of the term. Note 
that, due to the boundary conditions, each of the indi- 
vidual powers i,j,k>X. To lowest order in these new 
variables, the partial differential equation for H(u, v, w) 
is 



d 2 H 



d 2 H 



d 2 H 



uv- 



uw- 



vw- 







dudv dudw dvdw 

so that, when we put our lowest order terms in the equa- 
tion, the result is 



i+j+k—d 



(ij + ik + jk)cij^u l v : 'w k 



0. 



Since none of the exponents i,j or k are negative to en- 
sure that the solution has no additional poles - the very 
definition of H(u,v,w) ensures we have just enough to 
give a pre-classical solution - then we conclude that the 
coefficients Cj,j k = 0. Since this argument works regard- 
less of the value of the lowest order parameter d, the only 
solution we have to the differential equation with these 
boundary conditions is H(u,v,w) = 0. Therefore, the 
only pre-classical solution to the A = Bianchi I model 
is imi,m 2 ,m 3 = 0. When we include a non-zero cosmolog- 
ical constant, the equation to solve becomes more com- 
plicated, and there is no simple way to show that there 
is only the null solution. In this case, there is an incon- 
sistency between the coefficients of different orders in the 
series expansion of H(u,v,w), and so the only solution 
for the A/0 situation is also Cij,k = 0. 

These results are contrary to the expectation that the 
less constrained Bianchi I model would have at least as 
many, if not more, solutions than the isotropic case. One 
reason this does not occur may be in the chosen approach 
to quantization; a different method for quantizing the 
Hamiltonian constraint may allow more pre-classical so- 
lutions. For example, one can consider an alternate factor 
ordering when the constraint is quantized. In all the work 
done previously, such as the isotropic models, the triad 
operators are placed on the right 0, 0] . The problems 
seen in this present work may not carry over to another 
ordering, say a symmetric one. However, since we are 
considering a reduction of the full loop quantum gravity 
situation, one should use one choice of ordering consis- 
tently for all such symmetry reductions. If it turns out 
that our current choice is not the best for a specific model, 
this must be considered in the results of all other cases. 
This places a meaningful restriction on what methods 
can be used in the full theory of quantum geometry. 



On the other hand, it is possible that there are no tech- 
nical issues of this sort, and loop quantum cosmology is 
only feasible on certain types of space-times. Here, in 
the full Bianchi I model, there are no pre-classical solu- 
tions whatsoever. When we examine the Bianchi I LRS 
case 0, pre-classical states have been found, but they 
are rather limited. Near the singularity, the solution will 
assume a wave form that remains constant as the two pa- 
rameters m, n increase in value (here, m labels the two 
equal triad components, n the third). In fact, the only 
variation that occurs in the solution is in the n direc- 
tion; the sequence quickly assumes a constant value as 
one moves in the m direction of parameter space. Thus, 
for the two anisotropic cases considered so far where full 
quantum solutions are known, there are strong restric- 
tions on the possible states. It is possible that the inclu- 
sion of matter in the Bianchi I model may change this 
situation; a cautionary note would be that this does not 
occur in the closed isotropic model . It remains to be 
seen if these types of restrictions are solely for Bianchi 
I, or whether they occur more generically in anisotropic 
space-times. 
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